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We present a phase-resolved, optical, spectroscopic study of the eclipsing low- 
mass X-ray binary, EXO 0748-676 = UY Vol. The sensitivity of Gemini combined 



CN ' with our complete phase coverage makes for the most detailed blue spectroscopic 

00 ' 

j^^ ■ study of this source obtained during its extended twenty-four year period of activ- 

O ■ ity. We identify 12 optical emission lines and present trailed spectra, tomograms, 

P^ ' and the first modulation maps of this source in outburst. The strongest line 

CSl ■ emission originates downstream of the stream- impact point, and this component 

is quite variable from night-to-night. Underlying this is weaker, more stable 
axisymmetric emission from the accretion disk. We identify weak, sharp emis- 
r> ■ sion components moving in phase with the donor star, from which we measure 

_c^_' Keui = 329 ± 26kms~^. Combining all the available dynamical constraints on 

the motion of the donor star with our observed accretion disk velocities we favor 
a neutron star mass close to canonical (Mi ~ 1.5 M0) and a very low mass donor 
(M2 — 0.1 M0). We note that there is no evidence for CNO processing that is of- 
ten associated with undermassive donor stars, however. A main sequence donor 
would require both a neutron star more massive than 2M0 and substantially 
sub-Keplerian disk emission. 

Subject headings: binaries: X-ray, optical: stars, stars: individual: UY Vol 
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Introduction 



The structure and stability of compact stellar objects depends on the composition and 
equation of state (EOS) of their matter. Degenerate white dwarf matter consists of protons 
and electrons and can be described employing only special relativity, but in a neutron star, 
the nature of matter has changed and general relativistic effects are too large to be ignored, 
making the neutron star EOS a strong test of General Relativity. The EOS determines both 
the m aximum mass and the mass-radius relation of a neutron star (JLattimer fc Prakash 
2007h . 



The majority of neutron star masses measured in binary puls ars are consistent with the 
canonical value of 1.35 ± O.O4M0 (JThorsett &: Chakrabartylll999l ). however, low-mass X-ray 
binaries (LMXBs) have revealed a broader mass spread since accreting neutron stars grow 
beyond their birth masses. This i s consistent with higher masses deduced in millisecond 
pulsars (e.g., iDemorest et al.ll2010f ). which evolve from LMXBs. In order to discriminate 



between theoretical models and discern the appropriate EOS, it is important to accurately 
determine the mass distribution and the high-mass cutoff for neutron star stability. 

Because LMXBs are rich in emission lines, phase-resolved spectroscopy can reveal the 
mass and geometry of the system. Determining neutron star masses in LMXBs remains 
challenging, however, because most neutron star LMXBs are persistent X-ray sources in 
which the companion is never revealed directly. It has been shown in Sco X-1 and other 
sources that the motion of companion stars is visib le in fluorescent nitrogen and carbon lines 
(the 'Bowen blend' of N ill and C ill around 4640 A; iMcClintock et al.lll975). Several systems 
have shown sharp emission lines tracing the motion of the companion (ISteeghs &: Casares 



2OO2I : ICasares et al.ll2003l : iHynes et al.ll2003t ICasares et al.ll2006l : ICornelisse et al. 



2007a 



bijd). 



The ability to detect these narrow lines is dependent on the fraction of emission line light 
originating from the companion relative to that from the disk. 

Although knowledge of the radial velocity curve of the companion star is required for a 
complete determination of the system parameters, much can also be learned about the accre- 
tion disk by studying trailed spectra, Doppler tomograms, and modulation maps. Doppler 
tomography uses kinematic information from lines t o create a ve locity map of a system, even 
unraveling the complex information in blended lines ( lMarshll2005[ ) . Modulation mapping is an 
extension of D oppler tomography that permits mapping of time-dependent emission sources 
( lSteeghsll2003[ ). The analysis relaxes a fundamental constraint of Doppler tomography by 
allowing emission line flux to vary, then generates three tomograms simultaneously mapping 
the modulated and non-modulated emission. 



Among neutron stars, EXO 0748-676 (=UY Vol) has held particular promise for con- 
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straining the EOS of dense nuclear matter. UY Vol was first recognized by iParmar et al. 
( 119851 ) as a transient so urce, and then persisted in outburst for over twenty years before 



returning to quiesc ence f l Wolff. Rav. fc Wood II2008I : IWolff et al.ll2008l : iHvnes fc Jonesll2008 



Torres et al.ll2008l ). First identified with EXOSAT, UY Vol was later classified as a low- 
mass X-ray binary (LM XB) containing a neutron star because it shows Type I X-ray bursts 
( IGottwald et al.lll986l ). The source show s irregular X-ray dips and periodic X-ray echpses 
indicative of a 3.82 hour X-ray period (IWolff et al.l |2002| ). The X-ray eclipses require a 
high inclination, reducing a maior sour ce of uncertainty in the system parameter estimates 
( IParmar et al.lll986t iHynes et al.l 120061 ). The ephemeris of the source has also been studied 
extensively, revealing complex period changes (see, e.g.. IWolff et al.ll2002l]. though since th e 
source has recently transitioned into quiescence {Lx ~ 8.5 x lO^^erg/s; IWolff et al.ll2008l ). 



it is no longer possible to perform high-timing resolution analysis. There was considerable 
excitement associated with the puta tive detection of gra vitationally red-shifted absorption 
lines from the neutron star surface (ICottam et al.l 120021). Th is det ection could not be re- 
produced in a larger data set, however (ICottam et al.l l2008l ). and iLin et al.l ( l2010l ) have 
argued that these lines cannot originate from the neutron star surface in light of the now 
known 552 Hz spin frequency of the neutron star. Analyzing X-ray burst emission from the 



Ozell ( 120061 ) estimated a neutron star mass Mi > 1.8Mq and possibly as high as 



source. 

Ml = 2.10 ± O.28M0, although since this work used the gravitational redshift, this too may 

need reexamination. 



Wade et al.l (119851 ) identified the optical counterpart UY Vol, w hich brightened to V y 
17 mag during the outburst, and faded to R ~ 22 mag in 2009 (JHynes fc Jones! 120091 ). 



Pearson et al.l (120061 ) performed Doppler tomography on 15 optical and ultraviolet emission 
lines in UY Vol. Their study combined observations from HST, VLT, Magellan, CTIO 
4m, and lUE. Their data are consistent with a 1.35 Mq neutron star and main sequence 
companion, with a mass ratio q ~ 0.34, based on the gas s tream position shown in the 



Hell A4686 A Doppler tomogram. iMuiioz-Darias et al.l ( 20091 ) used intermediate resolution 
spectroscopy to perform Doppler tomography on narrow components in three He II emission 
lines. They measure i^em = 300 ± 10 km/s, and derive a neutron star mass Mi > I.5M0 
for the case of a main sequence companion. Modulation r aappin g was successfully applied 
to the quiescent counterpart of UY Vol by iBassa et al.l (l2009l ). Analyzing the Ha and 
He I A6678A lines, they directly identify narrow line emission from the counterpart and 
constrain Kcm > 400 km/s from their tomograms. They hence place a lower limit on the 
neutron star mass of Mi > 1.27Mq and estimate a m ass ratio 0.0 7 5 < q < 0.105 for the 
specific case of a 1.4 M© neutron star. Most recently, iRatti et al.l ( 12012| ) have performed 
further quiescent observations and obtain K^^ = 308.5 ±3.9 km/s from H/3 and H7 together 
with Fe II. Surprisingly, if the widths of their sharp components from the companion star 
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are due to rotational broadening, then vsini = 255 ± 22kms~^, and a neutron star mass in 
excess of 3.5 Mq is required. More likely is that the widths are non-rotational and signatures 
of material evaporating from the companion star in response to either X-ray irradiation or 
a pulsar wind, making this a 'Black- Widow-like' system. 

In this paper, we present optical spectroscopy of UY Vol obtained over four nights in 
2008, prior to the source's fade to quiescence. Our observations provide complete phase 
coverage over multiple orbits. We detect narrow emission components in the Bowen blend 



from the donor star, and find a value of ii'em consistent with the results of lMuiioz-Darias et al. 



(120091 ) . however, we do not confirm their observation of narrow Hell emission lines originating 
near the donor star. Our analysis represents the first attempt to deblend H l/He ll lines via 
Doppler mapping and also the first application of modulation mapping of this source in 
outburst. In Section 2, we describe our observations and data reduction. In Section 3 is our 
analysis and in Section 4, our discussion. 



2. Observations and Data Reduction 

We obtained intermediate dispersion, lon g slit spectroscopy with Gemini-South using the 



Gemini Multi-Object Spectrograph (GMOS: iHook et al.ll2004f ) on 2008 February 9-11 and 



March 14. The observations were designed to optimize spectral resolution for the observation 
of narrow emission components and to get phase coverage over multiple orbits. We obtained 
a total of 68 spectra each with 600s exposure, totaling 11.3 hours of time (> 3 binary orbits). 
With a 1.0 arcsec slit and B1200 grating, we obtain a resolution of 30 km/s at 4630 A. The 
spectra were reduced using the standard GMOS spectral reduction packages in IRAF|^. Over 
the four nights, we had an average seeing of 0.85 arcsec and a maximum seeing of 1.4 arcsec. 
Slit losses averaged 10%. Dome flats were taken nightly prior to and following each set of 
observations. The gemini.gmos data reduction package in IRAF generates normalized flat- 
fields, then apphes bias subtraction, fiat-fielding, and rough wavelength calibration using 
archival arcs. Wavelength calibrations are further refined using contemporaneous CuAr 
arc images taken during the daytime. Gemini gratings are calibrated daily, so the zero- 
points are well-known. Because Gemini has exceptional stability, the pixel-to-wavelength 
transformations are not affected by significant flexure and we did not require night-time 
arcs. The wavelength calibrated image is then sky subtracted, and spectra are extracted 



-'^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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using gsextract. Apertures were determined individually for each image, and variance 
we ighting is applied. To compute the phase for each spectrum, we use the ephemeris given 
by IWolff et al.l ( 20021 ). Our observations and phase coverage are listed in Tabled] 



We continuum fit each spectrum and generate both continuum normalized and con- 
tinuum subtracted spectra. We analyze both sets of spectra and verify that our choice of 
continuum treatment has negligible effect on our results. 



3. Analysis 

3.1. The Average Spectrum 

Combining 68 spectra, we create a continuum normalized, average spectrum shown in 
Figure [H We identify a total of 12 emission lines, sampling a range of ionization stages and 
excitations. These lines are listed in Table |2] with both their laboratory wavelength^! and 
measured line centers. 

Both Hi and Hell show double-peaked emission characteristic of orbital motion. To 
measure the observed line centers, we fit a double Gaussian to each line and determine the 
line center of each peak, then average those values together. By dividing our sample of 
68 spectra in half, we create two combined spectra, giving us two statistically independent 
measurements of the line center so that we can constrain our errors. The separation of each 
peak from the central wavelength for Hi and Hell is typically ~ 402.8 ± 0.5 km/s. 

We identify the Bowen blend in our sp ectrum. The Bowen b lend is a blend of Cm 
and Nlll lines ranging from 4638 - 4652 A (jMcClintock et al.l 119751 ). In typical LMXBs, a 
relatively equal proportion of C ill and N ill results in a b lended line observed a t 4640 A. Of- 
ten times, strong narrow emission features are observed ( ICornelisse et al.ll2008l ). In UY Vol, 
the Bowen line is centered at 4645 A, s uggesting a stronger co ntribution of Cm emission. 
GX 9+9 also shows this atypical ratio ( ICornelisse et al.ll2007bl ). We do not detect narrow 
line emission in our individual spectra. This likely means that the disk dominates the Bowen 
blend, but does not rule out weaker lines from the donor. 

Two He I lines A4922, A5016, a nd a blended O VI A5289/52 90 are also detected. This O VI 
feature was previously identified by lMuiioz-Darias et al.l ( l2009l ) as Fe ll A5284, but our trailed 
spectra (see below) clearly shows a movi ng line centered at A 5289. The identification of 
oxygen at this ionization is consistent with iPearson et al.l (120061 ) whose ultraviolet spectrum 



^The Atomic Line List v2.04 'http : //www . pa . uky . edu/~peter/atomic/ 
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showed V A1371 with very s imilar kinematics. Oth er LMXBs also show O vi line emission 
(see, e.g. ICasares et al.ll2003l : ISchmidtke et al.ll2000l ). We identify no other Fe lines in our 



spectrum of UY Vol, so O vi is a more plausible interpretation. 



3.2. Trailed Spectra 

Binning our 68 spectra by phase, we create trailed spectra for each observed line, and 
plot those in Figure |2l We observe clear modulation of the line flux with phase in all lines, 
and S-wave patterns emerge in several cases. The S- waves in the He ll trailed spectra have a 
hooked shape with a sharp dv/dcf) at < 0.5 and a smoother change for > 0.5. Emission is 
also relatively fainter at ~ 0.75. The Ovi and Bowen blend lines also show clear S-waves. 
The Ovi line shows higher emission velocity than H and He (710 km/s vs. 403 km/s), 
indicating either a location deeper inside the accret ion disk or high velo city disk overflow. 
Interestingly, in examining the LMXB 2A 1822-371. ICasares et allJ2003h find that the Ovi 
A3811 emission they observe has a radial velocity amplitude significantly larger than that of 
H/He components. 



Both He I and H l/He ll blends show characteristic absorption in the trailed spectra mak- 
ing it difficult to identify a clear S-wave pattern. In He I lines, absorption occurs consistently 
at 0.7 < (f) < 0.9. This is the same phase range at which X-ray dips are seen, suggesting the 
absorption may also be associated with the stream and its impact point. All lines showing 
S-waves have phasing consistent with emission originating in the accretion disk, either near 
the matter stream impact point or downstream of it. 



3.3. The Doppler Corrected Spectrum 



Mufioz-Darias et al.l (120091) suggested that sharp components could be seen in the Bowen 
line profile if it was Doppler-corrected to match a companion star radial velocity semi- 
amplitude of /Tem = 300 km/s, corresponding to a feature seen in their Hell Doppler to- 
mogram. We also find that sharp components can be recovered by Doppler correction at 
approximately this velocity, but have not been able to reproduce corresponding features in 
a tomogram. One should be cautious in interpreting results of Doppler correction in isola- 
tion, as there could be multiple spurious combinations of phase and radial velocity that will 
produce small apparent peaks out of the noise. 

We Doppler correct using the equation: 

1^(0) = 7 + ir^^ sin 27r(</.-</)o). (1) 
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To quantify the Doppler correction procedure more rigorously, we investigated the full param- 
eter space in phase (0-1 in steps of 0.02) and i^em (0-800 km/s in steps of 10 km/s). For each 
trial radial velocity curve we Doppler correct the Bowen profiles, and then cross-correlate 
them with a template consisting of one Gaussian matching the instrumental resolution for 
each N ill and Cm component. T h e indi vidual strengths of the Gaussians were set to match 
those listed by iMcClintock et al.l ( 1l975l ) with the overall strength of nitrogen and carbon 
emission adjusted to best match our data. Using a cross-correlation in this way should favor 
finding specific patterns of lines matching the Bowen blend wavelengths rather than random 
peaks in the data, and should also result in a single dominant peak for a good match, rather 
than the multiple peaks present in the corrected line profiles. 

After this procedure we find that the dominant feature in all of the cross- correlation 
functions is a broad peak corresponding to the broad emission associated with the disk. 
For phases near to zero and K^m near 300 km/s, a sharp peak is superposed on the broad 
maximum. To quantify the strength of this fit, we subtract a Gaussian fit to the broad 
maximum and then measure the height of the residual central peak. We show the strength 
of this peak as a function of the trial phases and K^m values in Figure O 

We see that the strongest feature is seen in phase with the motion of the donor star. The 
peak is at a phase of 0.98. We estimate i^em = 329 ±26 km/s with the 1 a uncertainty derived 
from a bootstrap Monte C arlo simulation of our data analysis procedure. The velocity found 
is consistent with that of iMuiioz-Darias et al.l (20091 ) and the va l ues de termined from the 
quiescent He I tomograms of iBassa et al.l (|2009[ ) and iRatti et al.l (120 121 ) . In Figure HI we 
show the Doppler corrected spectrum of UY Vol around the A4686 and 5412 He ll lines using 
the value ii'em = 329 km/s and 7 = 70 km/ s. We adopt the latter throughout this work for 
consistency with (JMuiioz-Darias et al.ll2009l ): we find a consistent value from our O vi line. 



While the Doppler corrected spectrum, shown in Figure H] shows reasonable centering 
around the line centers, we do not see the strong narrow line emission in Hell. We ob- 
serve at most a very smal l contribution to the He ll line from the companion star unlike 
Munoz-Darias et al.l ( 120091 ). ICornelisse et al.l ( 120081 ) summarize a survey of 10 LMXBs, all 
of which show He ll emission that is not localized. However, many of these sources show H I 
i n absorption. Both GX 9-1-9 and 2A 1822-371 sho w average spectra with Hi in absorption 
( ICasares et al.ll2003l : ICornelisse et al.ll2007bl . |2008|). All of our blended H/He hues show a 
dual-peaked shape when Doppler corrected, with one peak strongly centered near the Hell 
emission and decreased emission near H I. 



3.4. Tomograms 



We next utilize the continuum subtracted spectra to perform Doppler tomography, first 
rebinning the d ata to a cons tant velocity scale of 30 km/s/pix. We then create tomograms 
using doppler (JMarshI l2005l ) which utilizes a combined maximum entropy and chi-square 
minimization technique. In all cases, the entropy is measured relative to a Gaussian blur 
default, meaning entropy is insensitive to large scale variations. We generate tomograms 
extending ±1000 km/s around the line center. 

We use only the 49 spectra for which the source is out of eclipse. The maps, shown in 
Figure include model locations for the primary, secondary, matter stream, and LI point for 
a 1.5 Mq neutron star and an 0.1 M© donor (see Section [^^31) . The bulk of the line emission 
emanates from the accretion disk, especially downstream of the matter impact point, with no 
definitive contribution from the donor star. We were not able to reconstruct a tomogram of 
the Bowen blend. The weak narrow line emission and the volume of contributing components 
resulted in poor maps. 

In the top two rows of Figure [5], we show tomograms of our three H l/He ll lines. We 
used doppler to deblend the H/He contributions, and find that while the Hi contribution 
is distributed throughout the disk, the Hell contribution is concentrated just downstream 
of the matter impact point. The location of the deblended Hell emission is consistent with 
that seen in the isolated He ll lines. 



Using VLT data taken in 2003. |Pearson et al.l (120061 ) find He ll emission in A4686, 5412 is 
brightest s ignificantly clo s er to the matter stream and assumed location of the secondary. In 
both lines, iPearson et al.l (120061 ) observe concentrated emission with Vy > km/s. Our data 
very consiste ntly shows Vk ~ km/s. Our data do not show the matter stream emission 
observed by iMunoz-Darias et al.l ( l2009l ). Their data do show some emission features at 
negative Vy. Since their data were taken one month prior io ours, this precludes the existence 
of a concentrated source of He ll emission migrating along the streaming disk. 

Because our data set covers multiple orbits, we were able to group our spectra and 
characterize the changes in the tomogram on a timescale of days. In Figure O we show 
the trailed spectra and tomograms for three consecutive days for which we had sufficient 
orbital coverage. On February 9, the emission in the trailed spectra is notably brighter 
and the peaks broader. In the tomogram, this corresponds to emission concentrated just 
downstream of the matter impact point with a large spread of emission along the lower 
left tomogram. The next evening, the peaks in the trailed spectra are notably narrower 
and the corresponding tomogram emission is more compact. On February 10, the emission 
disperses toward the upper right and on February 11, the emission is most concentrated. 
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These variations in morphology also appear in the modulation maps (see below). Given how 
rapidly the centroid and morphology of He ll emission can change from epoch to epoch, it is 
unlikely that emission features in the tomogram associated with stream-disk impact r egion 



can be used to constrain the mass ratio reliably as attempted by iPearson et al.l (|2006[ ). 



Our He ll 4686 A maps do rather consistently define the accretion disk emission, and this 
potentially provides an alternative constraint on the system parameters. We show the same 
night-by-night tomograms in Figure [TJ with the gray-scale adjusted to maximize visibility 
of the disk on the right hand side of the tomogram. We will concentrate on this region as 
it should not be distorted by the changes seen in the stream-impact /bulge region that vary 
from night-to-night. Instead the kinematics of the right hand side of the maps are very 
stable and define a circular ring quite well. The projected disk velocity indicated by the 
data is 500 km/s, and is well bounded by rings at 400 and 600 km/s. We will discuss the 
implications of interpreting this as a Keplerian disk velocity in Section 14.31 

In the last row of Figure [5l we see that He I emits more strongly in a different region of 
the disk to Hell, separated by about 0.25 in phase. It is possible that the Hell gets ionized 
near the matter impact point and is able to later recombine and form He I downstream. 

The Ovi emission is fairly concentr ated, suggestin g temp erature conditions unique to 



O VI excitation at this point in the disk. iPearson et al.l (|2006[ ) observe O v emission in the 



UV. As with our O vi observations, their O V observations do not occur along the the ballistic 
stream or at Keplerian velocity , however, it is consi stent with a location downstream of the 



impact or the stream overflow (IPearson et al.l |2006|) . The high velocity of the O vi emission 



observed in the tomogram is expected from the fact that the S-wave seen in our O vi trailed 
spectrum had a higher amplitude than that of the H/He lines. 

The Doppler tomograms begin to form a consistent picture of the ionization structure 
of the disk. Hell is observed in the ionized disk bulge, predominantly on the left of the 
tomogram at negative Vx and low Vy- While we do observe some variance in the central 
location, we can say that we do not observe Hell emission from the donor star. We also 
know that the prominent Hell emission occurs downstream of the matter impact point. He I 
is absent from the ionized disk bulge. It occurs most prominently in the lower part of the 
tomogram at negative Vy, along the flow of the disk stream. As it occurs on the opposite 
side of the disk from the matter stream, it is probable that the temperature conditions in 
the disk are more amicable to the recombination of Hell. The Ovi emission, in contrast, 
may be tracing a higher ionization region inside the disk. 
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3.5. Modulation Maps 



Modulation mapping is an extension of Dopp ler tomograph y that maps sources varying 



harmonically as a function of the orbital period ( jSteeghsl |2003[ ) . Periodic absorption, disk 
flickering, and varying visibility of the emission region all contribute to features in the modu- 
lation map. Thus, we must go back to the trailed spectra for a more accurate interpretation. 
At present, the modulation mapping code does not account for eclipsing, so we use only 
49 of the 68 spectra for which the source was out of eclipse. Also, since the code does not 
deblend lines, we focus our study on isolated lines. Although modulation mapping works 
best with flux calibrated spectra, we did not have a comparison star on the slit for relative 
flux calibration. We assume that there is no significant continuum flux variation night to 
night and create modulation maps of all identified lines and plot two representative samples 
in Figures |8]|9l Modulation mapping considers a line source modulating harmonically with 
flux as a function of orbital phase, such that /(0) = lavg + I cos cos{2TT(f)) + /smsin(27r0). 
The modulation maps show the average, non-varying line flux {lavg) as well as the total 
modulated emission {Ism + Icos), and the individual sine and cosine emission maps. Because 
the modulated emission is often small relative to the average emission, the maps are shown 
on a fractional gray-scale relative to the average emission. Note that Isin and Icos can have 
positive or negative amplitudes and so their combination can model all phasings of the light 
curve. 

For the A4686 He ll line, the modulated emission contributes less than 6% of the average 
emission. Figure [8] shows maps of the He ll A4686 line. As with the Doppler tomogram, we 
observe the classic disk structure in emission, and concentrated emission just downstream 
of the stream impact point. Modulated emission is dominated by the sine term and appears 
most strongly downstream of the impact point. This corresponds to the region where the 
night-to-night tomograms differ most (see Fig. [6]). This modulated component would make 
our line emission appear brighter at ~ 0.25 and fainter at ~ 0.75. In the trailed 
spectrum, this is observed as a relative brightening in the S-wave near the peak radial 
velocity. A negative amplitude modulation appears in the cosine map, bordering the side of 
the disk farthest from impact. This would make our emission fainter at ~ and brighter at 
(p ~ 0.5. In the trailed spectrum, this appears as increased emission at ~ 0.5 in anti-phase 
with the S-wave amplitude. The modulating component may correspond to the irradiated 
inner edge of the stream bulge. 

The split structure of the He I A5015 emission is readily apparent in the modulation map 
shown in Figure [H Modulated emission here contributes up to 35% of the total flux. The sine 
term dominates the emission nearest the matter impact point. Very likely, this modulation 
does not represent varying visibility of an emission source, but rather the periodic absorption 
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which is observed in the trailed spectra at </> ~ 0.75. The non-modulated emission, lavg is 
strongest on the far side of the disk, where He ll is likely recombining. 

The presence of modulating emission contributes to our emerging picture of the disk. 
We certainly confirm the presence of ionized He ll emission at or near the ionized disk bulge. 
The modulating feature in the sine map could represent the formation of optically thick 
clouds above and below the plane of the disk. These clouds are formed when the matter 
stream impacts the disk, and they are illuminated by X-rays from the disk. Since they are 
optically thick to X-rays, we would only observe these clouds at ~ 0.25. These clouds may 
also be responsible for the periodic a bsorption in H e I and Hi and would be responsible for 
irregular X-ray dipping discovered by iParmar et al.l (Il986l ) . 



Discussion 



4.1. Accretion disk structure 



While Doppler tomography has long been used to interpret the accretion disk structure, 
we must be cautious not to over interpret certain aspects, as we have shown that prominent 
emission region morphology can vary between orbits . Ou r tomograms of Hell A4686 and 
5412 look different than those of iPearson et al.l (120061 ) and iMuiioz-Darias et al.l (120091 ). but 
this is likely reflective of the disk evolving between epochs. 



Pearson et al.l (120061 ) attempted to constrain the mass ratio of UY Vol by studying 



emission from the matter stream as observed in their tomograms for Hell A4686. Our 
observation of the same line does not show prominent emission from the matter impact 
point, but rather, some point on the disk downstream of the matter impact point. At best, 
our Hi maps show some emission coincident with the velocities expected for the matter 
stream, and He II is observed further downstream on the left side of the disk. The presence 
of He I emission in the lower part of the disk suggests the ionized He ll is recombining as it 
moves around the disk. 



Neither iPearson et al.l (120061 ) nor this work confir m the concentrated Hell e mission 
coincident with the location of the donor star reported by lMuiioz-Darias et al.l (120091 ). While 
the morphology of prominent disk emission is expected to change rapidly, emission from the 
donor star should be more consistent. It is possible that the relative brightness of narrow 
emission from the donor compared to the disk made it unobse rvable at other epochs. We 
note that the He ll emission associated with the donor star by iMufioz-Darias et al.l (120091 ) 
is part of a complex extending to negative Vx in their maps, and rou ghly parallelling the 
stream trajectory. It is then also possible the the Hell seen at \4 = by lMuiioz-Darias et al. 
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( 120091 ) originated from the beginning of the matter stream rather than the donor star itself, 
although it is unclear how He ll emission is produced in this region that should be shielded 
from direct irradiation by the disk rim. 

As we must already use caution when comparing tomograms of the same source from 
different epochs, so too we must cautiously approach ho w we compare th e tomo grams of 
UY Vol to other LMXBs. In examining the source GX 9+9, ICornelisse et al.l ( l2007bl ) observe 
prominent He II emission in the lower le f t qua drant of their tomograms, where UY Vol 
shows more prominent He I. iHynes et al.l ( 120011 ) find a similar emission structure for Hell 
in Doppler maps of XTE J2123-058. Both authors draw comparisons between these sources 
and SW Sex cataclysmic variables. The ionization structure of the emission source suggests 
a matter stream not c onnected with the accretion disk, that is ei ther overflowing the disk or 
being propelled away ( IHynes et al.ll200ll : ICornelisse et al.ll2007bl . and references therein). 



In UY Vol, the six He ll lines for which we create tomograms show emission most con- 
centrated at Vy ~ 0, but the Vx velocities occasionally suggest that while most emission is 
consistent with the Keplerian velocity of the outer edge of disk, some may be sub-Keplerian 
representing a similar overflow. 



Casares et al.l ( 120031 ) performed a study of the LMXB 2A 1822-371. Their map of He I 



A4471 places the He I absorber on the leading side of companion's Roche lobe or over the 
gas stream. It is not surprising that their He I maps are remarkably different than ours 
since they observe He I in absorption whereas UY V o l show s it in emission. Like UY Vol, 
2A 1822-371 shows oxygen emission. ICasares et al.l ( l2003l ) observe Ovi 3811 A clearly in 
the post-shock region between the matter stream and the disk. This is in a very different 
location than that observed for UY Vol. Our O vi emission emanates downstream of the 
matter impact point, at a high velocity, suggesting it occurs either deeper within the disk or 
in a region of high-velocity disk overflow. It is worth nothing that 2A 1822-371 is also at a 
higher inclination and likely higher mass accretion rate than UY Vol, thus the direct stream 
impact may be more excited. 

In summary, while UY Vol does not show the abundance of He I absorbers seen in 
2A 1822-371, it does seem to show similar ionized emission in the disk and possible high- 
velocity and sub-Keplerian disk overflows as seen in other LMXBs. Interestingly, in UY Vol, 
we observe He I emission as part of the accretion disk structure, which may give us some in- 
sight to the relative mass accretion rate and disk temperature as compared to other LMXBs. 
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4.2. Abundances 



We find that the C ill emission in the Doppler corrected Bowen blen d profile in Figure 
is of comparable stre ngth to Nlll. This can be compared to Sco X-1 (jSteeghs &: Casare; 



(jCornelisse et al.l l2007a[). anc 
4U 1735 —444. JCasares et al" 



GX 9+9 f Cornehsse et al. 



20021). 2A 1822-371 fICasares et al.M2003h. 4U 1636-536 fICasares et al.M 2006h. LMC X-2 



Aql X-1 (JCornelisse et al.l l2007cf ) where Nlll is dominant 



20061 ) where C ill and N ill are of comparable strength, and 



2007bl ) where C ill dominates. Thus the C lll:N ill ratio in UY Vol 



is among the higher values seen among neutron star LMXBs, suggesting that the carbon 
abundance is at least normal, if not enhanced. We should be wary of over-interpreting the 
ratio, however, as it does depend on physical conditions as well as a bundances. This i s 
illustrated by the diversity of Bowen blend structures seen in GX 339-4 (JHynes et al.ll2003l ). 



As discussed in the case of GX 9+9, the strong C ill rules out substantial CNO processing 
of the accreted material since this would result in a pronounced deficit of carbon relative to 
nitrogen. The effects of CNO processing on LMXB spectra are even more pronounced in 
the far-UV, where the normally dominant C iv line can become completely undetectable, as 
in the case of the black hole system XTE Jll 18+480 JHaswell et al.ll2002h . In UY Vol, on 



the o ther hand, the far-UV spectrum is typical with a very strong C iv line ( jPearson et al. 
2006h . 



4.3. System parameters 

There have now been four independent detections of features that appear to be asso- 
ciated with the companion star in UY Vol, spanning four different emission lines and both 
out burst and quiescence. In ou tburst we find i^em = 329 ± 26kms~^ in the Bowen blend, 
and iMunoz-Darias et al.l (120091 ) obtain K^m = 300 ± lOkms^^ from Hell. In quiescence 
Bassa et al.l ( l2009l ) obtain 345 ± 5 km s~^ in He I and 410 ± 5 km s^^ in Ha, while iRatti et al. 



( 2OI2I ) find 308.5+3.4 km s ^ from a weighted average of H/3 and H7. There is some spread in 



the deduced values, but there is also a striking congruence between them with emission from 
the companion star apparently spanning at least the range in velocity of 300-400 km s~^, 
with different average velocities seen at different times and in diffe rent emission lines. This 
provides more information than a single detection, as discussed by iBassa et al.l ( l2009l ). but 
also confuses interpretation as the 'K correction' from the observed Kem to the true K2 value 
( jMuiioz-Darias et al.ll2005l ) cannot be single valued. Part of the variation can be accounted 
for with ti me-dependent variat i ons in the disk shielding angle, a, which are included in the 
models of iMunoz-Darias et al.l (120051 ) but the iBassa et al.l ( l2009l ) measurement of different 



Kem values simultaneously in Ha and He I must indicate different patterns of emission across 
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the donor star. 

These limitations notwithstanding, we can attempt to constrain the masses of the bi- 
nary components u si ng th is in formation. In sona e respects we follow the discussion of 



Munoz-Darias et al.l ( l2009l ) and iBassa et al.l (120091 ). although we present the constraints 



more directly in M2 vs. Mi space in Figure [TOl The first constraint is that the K2 value 
must be greater than all measured -ft'em- The Ha detection is most constraining in this re- 
spect, requiring K2 > 405kms~^. We show this constraint in Figure [TUl as curve la which 
provides a lower limit on Mi as a function of M2. We can go beyond this limit if we assume 
that the emission originates in excitation by direct emission from near the compact object 
since the pole of the companion star itself (corresponding to Kem = K2) is shielded by the 
curve of the donor star's surface. A tighter constraint is that i^em is less than the velocity at 
the tangent point, where rays from the compact object are tangential to the stellar surface. 
We show this constraint as curve lb. 

While there is uncertainty in the actual emission pattern over the donor star, the center 
of light must be between the Li point and the tangent point, so we also require that i^Li be 
less than the smallest measured -ft'em, i-e. -f^Li ^ 310kms~^. This constraint is shown as 
curve 2 in Figure [TO] and this sets an upper limit on Mi as a function of M2. If Mi lies above 
this line then the mass ratio of the system is small and the donor is consequently too small 
to be consistent with the range of velocities seen by different groups at different times. 

These constraints together define a band in Mi — M2 space that becomes increasingly 
na rrow for sma l l mass ratios. This argument (using curves la and 2) is essentially that used 



by iBassa et al.l ( 120091 ) to constrain 0.075 < q < 0.105 (0.1 < M2 < 0.15) for the specific case 
of a 1.4 neutron star. This case is barely consistent with these constraints and requires an 
unusually extreme mass ratio and a very low-mass donor. Limiting emission to the tangent 
point rather than the pole (the band defined by curves lb and 2) virtually eliminates a 1.4 
neutron star from consideration. 

A 1.4 M0 neutron star is often assumed in LMXBs. There is some support for this 
i n the observed tight distributi on of masses of radio pulsars in binaries, of 1.35 ± 0.04 M0 



(iThorsett &: Chakrabartylll999l ). This tight distribution presumably reflects the distribution 
in birth masses, however, whereas neutron stars in LMXBs are accreting and so are likely to 
be of higher mass. Indeed, we do see rather larger masses where we can measure them in mill i- 



second pulsars which are believed to be the descendants of LMXBs (JDemorest et al.ll2010l ). 
Moving on to the companion star, the default assumption is often that it is a main-sequence 
star, or at least lies close to a main-sequence mass-radius relation. While under-massive 
mass donors are seen in some LMXBs, a main-sequence donor is still a reasonable starting 
point in the absence of other evidence for a more evolved status. Such other evidence 
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can sometimes be very strong, for example, in abundance anomalies indic ating that the 



dono r has experienced a prior period of CNO burning (e.g. XTE J1118+480; iHaswell et al. 



20021 ). In the case of UY Vol, however, as discussed in Section 14. 2^ it appears that the 



carbon abundance is at least normal, if not enhanced. If we follow the usual prescription 



for deducing the mass of a main-sequence donor in an interacting binary (JFrank et al.ll2002l ) 
we obtain M2 = 0.42 M©. This is based on very simplistic assumptions {R/Rq = M/Mq 
on the lower main-sequence), and more realistic mass radius relations produce somewhat 
lower don o r mas ses. By way of example, we show in Figure [TD] the mass-radius relation of 



Patterson! ( Il984j ) as curve 3. This favors masses around 0.37 M© for plausible neutron star 



masses. 

We can now draw one quite strong conclusion from the dynamical evidence in Figure [TUl 
if the donor star is on the main- sequence, then the neutron star is massive. If curve lb is 
obeyed (i.e. quiescent Ha emission is driven by direct irradiation from the neutron star or 
inner disk), then a main-sequence companion star requires a neutron star more massive 
than about 2M0. This case would actually agree quite well with the mass estimate of 



Ozell ( I2OO6I ) of 2.1 ± 0.28 Mq. Alternatively if we expect the neutron star to be closer to 
the canonical 1.4 Mq, then we require a significantly under- massive donor. In that case, 
we might expect to see evidence of CNO processing as in, for example, XTE J1118-I-480 



( IHaswell et al.ll2002l ). and we clearly do not. In this case, however, the low mass may be a 



consequence 



of ev aporation of an initially sub-Solar donor in a 'Black- Widow-like' scenario 



( JRatti et al.ll2012[ ) where no CNO processing has had timeto occur. 



We have one final dynamical constraint that has not been included yet. We observe 
the peak disk emission in Hell around a projected Keplerian disk velocity of VkepSini = 
500kms~^, and it is quite well bounded by a range of 400-600 km s~^. This turns out 
to be difficult to accommodate simultaneously with the dynamical constraints from the 
donor star. For a given pair of Mi,M2 values, the mass ratio and inclination are uniquely 
determined (via X-ray eclipses), and so the projected Keplerian velocit y from the disk rim 



can b e calculated. We assume the disk is tidally truncated at 0.9i?iobe ( IWhitehurst fc King 



19911 ). and refer to this velocity as ykep, tidal. This should be the minimum allowed velocity 
in a Keplerian disk. We show in Figure dU] a final set of curves (4a-d) corresponding to the 
loci Vkep,tidai sin i = 400, 500, 600, and 700kms~^ respectively. The 500kms~^ curve is 
not consistent with the dynamical constraints at all, and even 600kms~^ is only consistent 
with the lowest Mi and M2 values, and requires that the Hell emission be mildly sub- 
Keplerian. Accommodating a main-sequence donor would require that the He ll emission be 
sub-Keplerian by about 50%. The dynamical constraints from the disk velocities thus favor 
relatively small masses, but we cannot robustly rule out higher masses since all plausible 
solutions require some degree of sub-Keplerian motion. 
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We note that there is a caveat to this last point. ISomero et aLl (120121 ) have recently 
reported a similar analysis of 2A 1822-371, where they find Hell emission with disk-like 
morphology, but sub-Keplerian velocity. In the case of 2A 1822-371 the parameters are 
better constrained, so this conclusion cannot be avoided. That object also appears to be 
accreting at a substantially higher rate t hat UY Vol, however , and with significant mass loss 
( JBayless et al.ll2010l : lBurderi et al.ll2010l ). ISomero et al.l ( 12012| ) argue that their He ll emission 
originates in a disk wind above the orbital plane. It is then quite likely that 2A 1822-371 is 
an anomaly, and that UY Vol (with lower mass transfer rate) should exhibit more normal 
disk emission. 



5. Conclusions 



We have examined optical spectroscopy of the LMXB UY Vol spanning several orbits 
and observing nights. Doppler tomography shows that the strongest line emission, seen 
most characteristically in He ll, originates in the vicinity of the disk rim downstream of the 
stream-impact point. This dominant component is observed to be quite variable in intensity 
and in how far around the disk it extends. The variations are seen both within a night, 
and from night-to-night. This component is superimposed on weaker, less variable emission 
from the accretion disk. The phasing of the variations is consistent with the He ll emission 
originating on the irradiated inner side of a disk bulge region. The higher excitation O VI 
line is seen at a higher velocity than He ll, s uggesting some stream overfiow. It quite closely 
resembles the O v line seen in the far-UV by lPearson et al.l (120061 ) in its kinematic behavior. 



Careful examination of the Bowen blend reveals N ill and C ill components moving in 
phase with the companion star. These are best seen by Doppler correcting the line profiles 
using radial velocity curves and they ar e maximized for Kpjn = 329 =b 26km/s, consisten t with 
independen t meas urements of K^m by iMuiioz-Darias et al.l (120091 ) , iBassa et al.l ( 120091 ) , and 
Ratti et al.l (|2012l ). Combining these dynamical constraints with the inferred disk velocities 
from our Doppler tomograms, we favor a neutron star close to the canonical mass. Mi ~ 
1.5 Mq and a very low mass donor star, with M2 ~ 0.1 Mq, in spite of the lack of evidence 
for CNO processing. Even this solution requires that the peak disk-emission be mildly sub- 
Keplerian. Higher masses are disfavored but cannot be securely ruled out. A main-sequence 
donor star (M2 ~ 0.37 Mq) would require Mi > 2.0 Mq and peak disk emission at velocities 
30 % lower than expected at the edge of the disk. 
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Fig. 1. — The average, continuum normalized spectrum of UY Vol. We identify a variety 
of H/He lines with double-peak structures indicative of emission originating in the accretion 
disk. The He ll 4-9 falls on a chip gap. Measured line centers are listed in Table |2j 



Table 1. Log of CMOS Observations 



Date Start Time (MJD) End Time (MJD) Start Phase End Phase Orbits 



2008-02-09 


54505.193 


54505.370 


0.269 


0.398 


1.11 


2008-02-10 


54506.103 


54506.231 


0.981 


0.790 


0.81 


2008-02-11 


54507.171 


54507.337 


0.650 


0.728 


1.04 


2008-03-14 


54539.019 


54539.049 


0.576 


0.772 


0.19 
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Fig. 2. — Trailed spectra of ten prominent lines in UY Vol. The He ll lines show clear S-wave 
patterns. The Bowen blend shows a broad S-wave feature with no narrow emission lines. The 
Ovi line shows a clear S-wave with an amplitude -K'cm = 710km/s. Trailed spectra for He I 
and blended H l/He ll lines reveal phase-dependent absorption. Sine waves are over-plotted 
to guide the eye, with amplitude of 410 km/s and a phase offset of 0o = 0.3. 
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Fig. 3. — The strength of the residual cross-correlation peak after subtraction of the broad 
component of the correlation corresponding to the disk. A phase of or 1 corresponds to a 
radial velocity curves in phase with the motion of the donor star. 
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Fig. 4. — We create a Doppler corrected spectrum using Kcm = 329 km/s and 7 = 70 km/s. 
While we may tentatively identify narrow features in the Bowen blend, we do not observe 
strong narrow features in He II. Note that He II A4542 occurs on a chip gap, so our Doppler 
reconstruction of that line is incomplete. 



-24- 






u 
o 



500 



13 -500 

> 







o 
o 



500 




13 -500 

> 




OT 



o 

13 

> 



m 



o 

13 

> 



500 



-500 



500 



-500 



He II 
4100A 



He II 

4iggA 




He I 
4g22A 




H I 
4340A 





He II 
4338A 





He II 
4686A 




He I 
5015A 




H I 
4861 A 



He II 
4859A 



He II 
5412A 




VI 
5289A 



-500 500 
Velocity (km/s) 



-500 500 
Velocity (km/s) 



-500 500 
Velocity (km/s) 



Fig. 5. — Doppler tomograms of identified lines. The components for the Hi and Hell hnes 
are deblended and shown separately. 
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Fig. 6. — Night-by- night tomograms of Hell 4686. Note, the tomograms for 2/10 and 2/11 
are shown at a fractional gray-scale relative to the 2/9 for clarity. They only achieve 50% of 
the signal. This is apparent by the relative peaks in the trailed spectra. 
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Fig. 7. — Night-by-night tomograms of He ll 4686 as for Figure [6l The gray-scale has now 
been adjusted individually to maximize visibility of the disk on the right hand side of the 
tomogram. The plotted circles correspond to disk velocities of 400, 500, and 600kms~^. 
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Fig. 8. — Modulation maps of the He ll A4686 line. Starting at the upper left and going 
clockwise, the panels are (1) the average non-modulating emission, (2) the total modulated 
emission {Igin + hos), (3) the Isin component of the modulated emission, and (4) the Icos 
component of the modulated emission. The modulation maps are shown at a fractional 
gray-scale of ±6% of the average, non-modulating emission. 
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Fig. 9. — Modulation maps of the He I A 5015 line. The panels are as described inO The 
fractional gray-scale employed here is ±35% of the average non-modulating emission. 
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Fig. 10. — Dynamical constraints on the mass of the neutron star (Mi) and the companion 
star (M2). The horizontal solid line corresponds to a canonical 1.4 (ti neutron star. Th e solid 
lines la and lb are derived from the quiescent Ha measurement of iBassa et al.l (l2009l ). The 
lower line, la is a lower boundary of allowed parameter space assuming only that K2 > -ft'em- 
The higher line, lb requires additionally that the Ha emission originates in regions directly 
ill uminated by the neutron st ar. The dashed line, 2, is obtained from the Hell measurement 
of iMuiioz-Darias et al.l (120091 ) . assuming that i^Li < -f^em- The dot-dashed line 3 represents 
main-sequence donors that fill the companion star Roche lobe. Finally, the dotted lines are 
constraints imposed by our measurements of the disk velocity. Curves 4a-d respectively 
correspond to outer disk velocities of 400, 500, 600, and 700 km/s. 
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Table 2. IDENTIFIED LINES 



Line 


^air 


Ae(A) 




Hi 2-6 
Hen 4-12 


4101.734 
4100.041 


4101.6 ± 0.6 


Hen 4-11 


4199.832 


4201. ± 2. 


Hi 2-5 
Hen 4-10 


4340.464 
4338.671 


4340.8 ± 0.9 


He II 4-9 


4541.591 


- 



Cm 4647.4 4645.0 ± 0.8 

4650.3 

4651.5 
Nin 4634.2 

4640.6 

4641.9 



He II 3-4 


4685.71 


4686.3 ± 0.8 


Hi 2-4 
He II 4-8 


4861.352 
4859.20 


4861.4 ± 0.7 


He I 


4921.93 


4922. ± 1 


He I 


5015.67 


5015.3 ± 0.8 


C/N/0 blend 


- 


5133 


1 — 1 1— 1 

> > 

O O 


5289. 
5290. 


5289. ± 1. 


He II 4-7 


5411.53 


5412. ± 1. 



Note. — Table of identified lines. The second column is the laboratory predicted central 
wavelength in air. The third column is the best estimate of the centroid as measured from 
the average spectra. Errors are estimated by creating two statistically independent spectra 
each summing over 34 of our nights. The A4541 He ll line falls on a chip gap, so no center is 
measured. The feature visible at 5133A is likely a blend of C/N/0. 



